The Josephson critical current in a long mesoscopic S-N-S junction 



P. Dubos, H. Courtois, B. Pannetier 

Centre de Recherches sur les Tres Basses Temperatures-C.N.R.S. associated to Universite Joseph Fourier, 25 Ave. des 

Martyrs, 38042 Grenoble, France 

F. K. Wilhelm,i'2 A. D. Zaikin^ and G. Schoni 
^ Institut fiir Theoretische Festkorperphysik, Universitdt Karlsruhe, 76128 Karlsruhe, Germany 
^Quantum Transport Group, Department of Applied Physics and DIMES, TU Delft, 2600 GA Delft, The Netherlands 
^ Forschungszentrum Karlsruhe, Institut fiir Nanotechnologie, 76021 Karlsruhe, Germany 

(February 1, 2008) 

We carry out an extensive experimental and theoretical study of the Josephson effect in S-N-S 
junctions made of a diffusive normal metal (N) embedded between two superconducting electrodes 
(S). Our experiments are performed on Nb-Cu-Nb junctions with highly-transparent interfaces. We 
give the predictions of the quasiclassical theory in various regimes on a precise and quantitative 
level. We describe the crossover between the short and the long junction regimes and provide 
the temperature dependence of the critical current using dimensionless units eRj^Ic/ec and ksT/tc 
where ec is the Thouless energy. Experimental and theoretical results are in excellent quantitative 
agreement. 
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The Josephson effect is well known to exist in weak 
links connecting two superconducting electrodes S, e.g. 
a tunnel barrier I, a short constriction C or a normal 
metal N (S-I-S, S-C-S and S-N-S junctions). This effect 
manifests itself in a non-dissipative DC-current flowing 
through the Josephson junction at zero voltage. At weak 
coupling, e.g. in the S-I-S case, the Josephson current 
can be expressed as Is = Ic sin ip where (p is the phase 
difference between the two superconducting condensates 
and the maximum supercurrent Ic is called the critical 
current. 

The Josephson effect in S-N-S junctions has been stud- 
ied in a V|ariety of configurations. The early experiments 
of Clarktfl and Shepherda were performed in Pb-Cu-Pb 
sandwiches. In these experiments and in the pioneer- 
ing calculations by de Gennes,u it was already realized 
that the presence of a supercurrent in such structures is 
due to the proximity effect. It can be understood as the 
generation of superconducting correlations in a normal 
metal connected to a superconductor, mediated by phase- 
coherent Andreev reflections at the S-N interface. The 
critical current Jc is limited by the "bottleneck" in the 
center of the N-layer, where the pair amplitude is expo- 
nentially small, Ic oc e^^/^'^. Here Lt = y^hD/2nkBT is 
the characteristic thermal length in the diffusive limit and 
L is the length of i-he junction. These calculations, as well 
as those by Fink,cl analyzed the temperature dependence 
of Ic within the Ginzburg-Landau theory in the vicinity 
of the superconducting critical temperature Tc. pLater, 
the critical current of diffusive S-N-|S microbridgesoQ was 
successfully described by LikhacevEl with the aid of the 
quasiclassical Usadel equations.l In this work, the em- 
phasis was put on the high temperature regime where 
the superconducting order parameter is smaller than the 
thermal energy A <C kBT. A more general study of the 
Josephson effect in diffusive S-N-S junctions was made in 



Ref. |. 

More recently, experimental data on long Joseph- 
son junctiongla showed a surprising temperature- 
dependence, which turned out to be in a strong disagree- 
ment with the early theory by de Gennes. These data 
have been discussed by some of usllil within the quasi- 
classical anproach which we will also use in-Lhe present 
work. Finkllj attempted to analyze the datatj by means 
of an extrapolation of the Ginzburg-Landau theory to 
low temperatures. 

The proximity effect in mesoscopic hybrid structures 
consisting of normal and superconducting metals at^ 
tracted a growing interest during the recent years .Ej 
Here we will consider mesoscopic diffusive S-N-S junc- 
tions where the sample length is much larger than the 
elastic mean free path 1^ but smaller than the dephas- 
ing length : < L < L^. In N-S junctions and 
Andreev interferometers, we can identify - both theo- 
retically and experimen t all) ; - the natural energy scale 
for the proximity effect .oEJ It is given by the Thouless 
energy Cc = TiD/L"^. Here D = wf4/3 is the diffusion 
constant of the N-metal, vy is the Fermi velocity. In con- 
trast to the energy gap A which is set by the interactions 
in the superconducting electrodes, the energy scale Cc is a 
single-electron quantity : Cc/?* is merely the diffusion rate 
across the sample for a single electron. This energy scale 
remains important in non-equilibrium situations, e.g. if 
one drives the supercurrent across a S-N-S junotkm^by 
the injection of a control current in the N-metalJl3il3 

The main purpose of the present paper is to carry out 
a detailed experimental investigation of the equilibrium 
supercurrent in relatively long diffusive S-N-S junctions 
with highly transparent N-S interfaces as well as a quan- 
titative comparison of our data to the theoretical pre- 
dictions. Here, a long junction means that the junction 
length L is much bigger than ^JfiD / IS.. This is equivalent 
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to A ^ ec. In order to perform this comparison at all rel- 
evant temperatures, we complete the previous studies by 
providing a rigorous expression for the Josephson critical 
current at T —^ which was not properly evaluated be- 
fore. Our experimental results arc in excellent agreement 
with the thccjccjtical predictions. 

As beforejBO our theoretical approach is based on the 
quasiclassical Green's functions in imaginary time. The 
proximity effect is described by a finite pair amplitude 
F in the N- metal (see [|l9) and references therein). We 
will assume N-S interfaces to be fully transparent and 
neglect the suppression of the pair potential A in the S 
electrodes near the N-S interface. This is appropriate at 
T <C Tc or if the reservoirs are very massive as compared 
to the normal metal. Within those bounds, our calcula- 
tion does not contain further approximations and is e.g. 
valid at arbitrary temperature and sample size. We will 
now proceed by discussing certain limits. 

In the high temperature regime k^T 3> Ec (or, equiva- 
lently, L ^ ^t), the solution is well known. In this case 
the mutual influence of the two superconducting elec- 
trodes can be neglected and the Usadel equations can be 
linearized in the N-metal,|-except in the vicinity of the 
N-S interfaces. One finds 
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where _Rn is the N-metal resistance, 
is the Matsubara frequency, fin — y/A^ 
y/KDj2uj^. If T is close to the critical temperature of 
S, the gap is small as compared to the thermal energy : 
A <C fceT. In this lirwt, Eq. (|l|) coincides with the result 
derived by Likharev.Q 
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FIG. 1. Calculated dependence of the zero temperature 
cRnIc product in units of A as a function of the ratio e^/ A. 
Ic is the Josephson critical current, i?N the normal state resis- 
tance, £c is the Thouless energy and A is the superconducting 
gap of S. The long junction regime is on the left part of the 
graph where tc < A, the short junction regime is on the right 
part where tc > A. The 4ftshed line corresponds to the Ku- 
Uk-Omel'yanchuk formulaEa at T = 0. 



At lower temperatures fceT" ^ Ec evaluation of Ic pr 
volves solutions of the Usadel equation sJ. all energiesE3. 
In order to determine the precise valuec^l of the critical 
current, we performed a numerical solution of the Usadel 
equations for the whole range of Matsubara frequencies. 
In the long junction limit (A 3> Cc), the zero-temperature 
eR^Ic is found to be proportional to ec : 



ei?N/c(T = 0) = 10.82ec 



(2) 



In this case, the current phase relation is slightly different 
from a sine and the supercurrent maximum ocr.iirs at 
(/? = 1.277r/2.As compared to previous estimates, ElHi3 the 
exact numerical prefactor in this formula turns out to 
be unexpectedly high. This observation is crucial for a 
quantitative comparison between theory and experiment 
not only in the case of copventional junctions but also 
for high-TcpSj-N-S junctionscJ or devices involving carbon 
nanotubesEj 

Let us now turn to the short junction regime A ^ Cc, 
i.e. to the case of dirty S-C-S weak links described in Ref. 

Our numerical results reproduces quantitatively 
the behaviours of both the current-phase relation and the 
zero-temperatur e cri tical current : eR^Ic ~ 1.3267rA/2 
at iy9 = 1.257r/2.E3E3 This results confirms the precision 
of our calculation in describing both the long junction 
and the short junction regimes. Our numerical results 
for Ic [T = 0) as a function of the Thouless energy ec are 
presented in Fig. |^. It confirms that it is the minimum 
of the gap A and the Thouless energy Cc which limits 
the critical current in diffusive S-N-S junctions. At tc — 
A, the critical current value remains close to the short 
junction case. 
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FIG. 2. Calculated temperature dependence of the bRnIc 
product. The different curves correspond to various values 
of the ratio A/tc in the long juction regime. The curve for 
A/ec ^ oo is universal in the sense it does not depend on A. 
Note that fcBT/ec = L^/27rL|. 

In the following, we will focus on long junctions A > ec- 
Fig. H shows the temperature dependence of the eRj^Ic 
product for various values of the superconducting gap 
in the long junction regime. Both axis are given in 
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units of the Thouless energy. The low-temperature part 
(fceT < 5ec) conies from a numerical solution of the 
Usadel equation while the high-temperature part comes 
from Eq. 0. From this figure, we can see that the char- 
acteristic decay temperature for the critical current is a 
few times the Thouless temperature ec/k-Q. As soon as 
k^T > 5ec, the sum in Eq. (|l|) can be reduced to the first 
frequency term within a 3% underestimation. This term 
corresponds to ujq = -Kk^T and L^^ — Lt- Adding the 
second term in the summation decreases the error below 
0.1% in the same temperature range. 

The universal curve of Fig. I for A/ec oo is valid 
only in the case of a very long junction with A/ec ^ 100. 
It appears as if A is to be compared to the quantity 
eR^Ic{T = 0) ~ lOcc in the long junction limit. In the 
limit of infinite A/ec, Eq. ^ simplifies to 
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From Eq. |3[ one can get the temperature dependence of 
the critical current : 1^ oc T^^''^exp{—L/Lj^). It has been 
demonstrated in Ref. [ ^ that within a limited temper- 
ature interval this expression is numerically very close to 
a simple exponential dependence Ic oc exp(— L/Lt) p4ft 
Lt 1/T, as one would expect in a ballistic limit.BLj 
From Fig. |, the quasi-exponential temperature depen- 
dence of the critical current is indeed striking. This 
was the central result of Ref. [ |T^, but was not un- 
derstood at that time. This coincidjCnce is purely acci- 
dental and has no special meaningO In the low temper- 
ature limit, the numerical solution can be approximated 
by ei?N/c/ec = a{l - 6e-°'='=/3-2feBT-)^ j,^^ coefficients a 
and b are 10.82 and 1.30 respectively in the long junction 
limit A/ec oo. 

S-N-S junctions are intrinsically shunted and have neg- 
ligible internal capacitance, so they are strongly over- 
damped. Their current-voltage characteristics are hence 
intrinsically non-hysteretic. The transition from a super- 
current to a voltage state happens at the ^critical cur- 
rent, but is rounded by fitute tcmperaturctZl We fabri- 
cated Nb-Cu-Nb junctionsES with a large conductance so 
that thermal fluctuations remain small compared to the 
Josephson energy : ksT <C hIc{T)/2e even at high tem- 
perature near the critical temperature of Nb. This makes 
a well-defined critical current up to the critical tempera- 
ture of Nb. Effects of environmental fluctuations known 
from mesoscopic tunnel junctions,t3 which are intrinsi- 
cally underdamped, are absent. _ 

We benefited from a trilayer stencil mask technologytB 
making use of a thermostable resist that does not outgas 
during Nb evaporation. Thus we were able to routinely 
obtain a superconducting critical temperature as high as 
S.IK for the Nb electrodes. We performed successive 
shadow evaporations of Cu and Nb at different angles 
through the silicon stencil layer in an ultra high vacuum 
chamber, followed by a lift-off. Fig. ^ shows a typical 
sample. We studied a single sample (a) plus five different 



samples evaporated on the same substrate (b, c, d, e 
and f). Table 1 lists the main physical parameters for 
these samples. The Cu metallic strips are 600 rtm wide 
and 100 nm thick. The Nb superconducting electrodes 
are 800 nm wide and 200 nm thick, except for sample a 
where it is 100 nm. The length L of the metallic island 
was varied between 700 and 1000 nm, corresponding to 
a separation length c^Nb between Nb electrodes varying 
between 370 and 700 nm. For all samples, the calculated 
Thouless energy fiD/L? is therefore significantly smaller 
than the gap A. 




FIG. 3. Oblique micrograph of a typical S-N-S junction 
made of a Cu wire embedded between two Nb electrodes. The 
doubling of every structure due to the shadow evaporation 
is visible. The Nb electrodes cover the Cu strip over about 
150 nm. 

The normal-state resistance i?N cannot be directly 
measured at temperature above Tc since the resistance of 
the Nb electrodes is measured in series. We found that 
the finite-bias resistance {eV ~ ec) varied by about 10% 
between 2 K and 8 K due to the proximity effect on the 
conductance. We took for the normal-state resistance 
Rn the resistance at T = 6 for a better agreement 
with the theory. It is a relatively high temperature since 
k^T > 15ec for every sample. Using L for the Cu length, 
we obtain a Cu resistivity p = 1.1 . 10~^ fJ.m for samples 
b to f and p = 1.5 . 10~^ fi.m for sample a. 

We measured the critical current of samples a to f at 
temperatures down to 300 mK. Our procedure consists 
in sweeping the bias current while measuring the dif- 
ferential resistance dV/dl. We define the experimental 
critical current as the current where the differential re- 
sistance reaches With this criteria, the experi- 
mental uncertainty is estimated below 0.5% at T ^2K, 
5% around T = 4is: and 100% at 7K. Fig. | shows 
the data for 3 samples. The measured eR^Ic/cc plotted 
as a function of the reduced temperature k^T / show 
a large decrease over more than two decades. For each 
sample, we fitted the data to the theoretical prediction 
with only one fitting parameter, the Thouless energy. 
The zero-temperature superconducting gap A was calcu- 
lated from the measured critical temperature of Nb using 
: A = 3.8fcBT'cEil This gives 1.3 mel/ for all samples ex- 
cept sample a for which A = 1 meV . We used both a 
fixed gap equal to the zero-temperature value and a gap 
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A(T) following the BCS temperature dependence, but 
with a slightly reduced critical temperature Tc = 7.5 K. 
At high temperature, it appeared necessary to take into 
account the temperature dependence of the gap. In this 
case, the agreement between theory and experiment is 
excellent. The fit is very sensitive to the chosen value of 
the Thouless energy. We would like to stress that for each 
sample the horizontal and vertical axis are normalized to 
the same Thouless energy ec- This fitted values are found 
to be very close to the Thouless energies calculated from 
the full length L of the Cu strip, see Table |. 
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FIG. 4. Temperature dependence of the measured bRnIc 
product of samples a, e and f together with the theoretical 
fits assuming a temperature-independent gap (full line) and a 
gap following a BCS temperature dependence with Tc = 7.5 K 
(dashed line). The only adjustment parameter is the Thou- 
less energy Ec of each sample. For description of the sample 
parameters, see Table 1. 

In Fig. |[ the critical current of sample f shows the 
onset of the saturation regime. At T = 30QmK the 
adjusted critical current eR^Ic reaches up to 8.2 ec- This 
number is close to the theoretical value 8.79ec for sample 
f at r = 0. This result discards an interpretation of our 
data within the Ginzburg-Landau theory of Ref. ^ which 
predicts a maximum eR^Id of about 1. 

In Ref. |l^, an array of S-N-S junctions was made of 
a long N-metal wire periodically in contact with a series 
of superconducting islands. A good fit between the data 
and the theory was shown in Ref. [ |Tl|, but with the 
introduction of a strong reduction of the effective area. 
This may be attributed to the periodic and lateral char- 
acters of this type of samples. 

Our calculation assumes perfectly transmitting inter- 



faces with zero boundary resistance. In_feLCt, it is suf- 
ficient that the barrier-equivalent lengthEj Lt = Ic/t is 
much smaller than the sample length. As an example, 
this condition means an interface transparency t > 0.1 
for sample b. In the case of Nb-Cu-Nb samples fabri- 
cated through a two-lithography-step process including 
Ar-etching,t3 we found a critical current with a reduced 
magnitude, presumably due to a slightly degraded inter- 
face. The critical current in S-N-S junctions with par- 
tially transparent interfaces was discussed in Ref. [ p4[ . 
The predicted behavior features a different temperature 
dependence for the critical current. Nevertheless, the 
temperature dependence remained consistent with the 
theory assuming a perfect interface. Only a reduction 
prefactor had to be introduced. This observation could 
hint at the fact, that interface barriers are very inho- 
mogenous and the current is carried through a few highly 
conducting pinholes. 

In summary, we discussed the Josephson critical cur- 
rent of diffusive S-N-S junctions. This study provides a 
simple and reliable formulation that enables the practi- 
cal determination of the equilibrium critical current. We 
studied the critical current of a set of samples with dif- 
ferent junction lengths and found an excellent agreement 
between our data and the predictions of the quasiclassical 
theory. 
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a 


1000 


600 


600 


0.260 


200 


13 


14.3 


70 


8.91 


b 


1010 


680 


580 


0.173 


300 


20 


18.6 


70 


8.99 


c 


910 


570 


590 


0.179 


260 


22 


21.7 


60 


8.83 


d 


800 


470 


580 


0.183 


230 


25 


25.4 


51 


8.64 


e 


800 


476 


590 


0.169 


250 


26 


26.1 


50 


8.62 


f 


710 


370 


580 


0.152 


250 


34 


33.5 


39 


8.32 



TABLE I. Parameters of the different samples studied. L is the full length of Cu strip while is the Nb electrodes 
separation, w is the Cu strip width. The Thouless energy tc is derived from the fit of the experimental data to the theoretical 
prediction (see Fig. 0). 
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